arXiv:1502.00282v2 [cond-mat.dis-nn] 28Jun2015 


Finite-Temperature Fluid-Insulator Transition of Strongly Interacting ID Disordered 

Bosons 


V.P. Michal 1 , I.L. Aleiner 2 , B.L. Altshuler 2 ’ 5 , and G.V. Shlyapnikov 1 ’ 3,4,5 
1 Laboratoire de Physique Theorique et Modeles Statistiques, Universite Paris Sud, CNRS, 91405 Orsay, France 
2 Physics Department, Columbia University, 538 West 120th Street, New York, New York 10027, USA 
3 Van der Waals-Zeeman Institute, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands 
4 Russian Quantum Center, Novaya street 100, Skolkovo, Moscow region 143025, Russia 
5 Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China 

(Dated: June 30, 2015) 

We consider the many-body localization-delocalization transition for strongly interacting one¬ 
dimensional disordered bosons and construct the full picture of finite temperature behavior of this 
system. This picture shows two insulator-fluid transitions at any finite temperature when varying 
the interaction strength. At weak interactions an increase in the interaction strength leads to 
insulator—>fluid transition, and for large interactions one has a reentrance to the insulator regime. 


In spite of intensive studies during several decades, An¬ 
derson localization of quantum particles in disorder [l] 
remains one of the most active directions of research in 
condensed matter physics 0. A subtle question is how 
interactions between the particles affect the localization. 
It was first raised for electrons in solids SI and is now 
becoming crucial for ultracold neutral atoms in random 
potentials SI- After the first experiments on the obser¬ 
vation of Anderson localization in expanding dilu te q uasi- 


one-dimensional (ID) clouds of bosonic atoms [9|,ll0|, the 
research on ultracold quantum gases in disorder rapidly 
grows. 

Interacting quantum particles can undergo many- 
body localization-delocalization transition (MBLDT) - 
the transition from insulator to fluid state | . MBLDT in 
a system of disordered weakly interacting ID bosons has 
been discussed both at zero EHH and finite tempera¬ 
tures |. In the latter case MBLDT manifests itself as a 
non-conventional insulator-normal fluid phase transition, 
the transport properties being singular at the transition 
point: in the fluid phase the mass transport is possible, 
whereas in the insulator phase it is completely blocked 
although the temperature T is finite. The fluid-insulator 
transition for strong interactions at T = 0 has been dis¬ 
cussed by Giamarchi and Schulz 14], who predicted that 
for the Luttinger liquid parameter K < 3/2 even an ar¬ 
bitrary weak short-range disorder leads to localization. 

Up to now the finite temperature behavior of ID disor¬ 
dered bosons was well understood only for weak interac¬ 
tions. In this Letter we extend this understanding to the 
general case of strong and moderate interactions. At any 
finite temperature, we show the presence of two insulator- 
fluid transitions: the insulator—>-fluid transition when the 
interaction strength increases from zero to a certain crit¬ 
ical value, and then a reentrance to the insulator phase 
at sufficiently strong interactions. 

The physical picture can be interpreted as follows. Lo¬ 
calization of all single-particle quantum eigenstates in ID 
by an arbitrary weak disorder 15. M implies the insulat¬ 
ing phase in the absence of the interaction between the 


bosons. In Ref. [/] it was demonstrated that arbitrary 
weak interactions are unable to destroy the insulator: the 
boson density gets fragmented into lakes with irrelevant 
tunneling between them. The tunneling becomes rele¬ 
vant and drives the system into a fluid state at a critical 
interaction strength determined by the disorder. 

On the other hand, it is well known that in the absence 
of disorder bosons with an infinitely strong repulsion are 
equivalent to free fermions 0 and, hence, they will be 
localized by an arbitrary weak disorder. Accordingly, 
there should be a second critical interaction strength 
above which one should expect an insulating state. In or¬ 
der to describe this transition using the ID boson-fermion 
duality one should determine an effective interaction be¬ 
tween the fermions when the boson-boson interaction is 
strong but finite. It is also important to account for the 
renormalization of the disorder by the fermion-fermion 
interaction even if the latter is fairly weak (see below). 

The system of N one-dimensional bosonic atoms repul¬ 
sively interacting with each other via a short-range po¬ 
tential is well described by the Lieb-Liniger model jl9|. 
The Hamiltonian reads: 

f.2 n 

Hb = -^J2 d x j + 9 X 8(xj-x k ), (1) 

j =1 1 <j<k<N 


where g > 0 is the coupling constant, and m the atom 
mass. The interaction strength is characterized by the 
dimensionless parameter 


_ mg 
7 h?n 


( 2 ) 


with n = N/L being the average density of bosons, and 
L the system length. In the regime of strong interactions 
we have 7 1, and in the Tonks-Girardeau limit of 

an infinitly strong repulsion (7 —> 00) the system maps 
onto the ID gas of free fermions 47]. In the general 
case of an arbitrary interaction strength, the Lieb-Liniger 
model m maps onto spinless (spin-polarized) fermions 
interacting with each other via an odd-wave momentum- 
dependent attractive interaction )[2ol - l23j . More precisely, 
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the fermions are governed by the Hamiltonian 


h 2 


N 




■ Vf, 


( 3 ) 


3 =1 


with the interaction operator [ 21 ] 
ft 4 


Vf = 


m 2 g 


Y (d Xj -d Xk )6{x j -x l: ){d Xj -d Xk ), (4) 


1 <j<k<N 


and the eigenfunctions of Hp 0 coincide with the 
bosonic eigenfunctions of Hb 0 when the coordinates Xj 
are ordered: Xj < Xj+ ±. Strongly repulsive bosons map 
onto weakly attractive spinless fermions with the Fermi 
momentum kp = 7 rn and Fermi energy Ep = h 2 k 2 F /2m. 
In the degenerate regime the fermionic dimensionless cou¬ 
pling constant can be estimated as the ratio of a typical 
interaction energy 0 per particle to the Fermi energy: 


A ~ — 


h 4 k 3 p/ 


m 2 g 


Ep 


-i«l. 

7 


( 5 ) 


Having in mind the comparison between strongly and 
weakly interacting bosons we use below the temperature 
of quantum degeneracy 

7ft = 


2 m 

which coincides with Ep of effective fermions in the case 
of strongly interacting bosons. 

Without loss of generality, we can represent the static 
disorder by a Gaussian random potential U (x) with zero 
mean, variance Uq, and correlation length cr: 

(U(x)) = 0; (U(x)U(x')) = Ulf{\x - x’\/a), 

where f(z) —>■ 1 for z —>• 0 , f{z) —>■ 0 for 0 —>• 00 , and 
the symbol (...) denotes averaging over the realizations of 
the disorder. Below we assume that Uq <C ftf/mcr 2 . The 
energy scale brought by the disorder is 25 h27I] 


£* = 


mcr 2 17g 

H 2 


1/3 


<t/n 


Derivation of this equation, which was made in [tJ , is de¬ 
scribed in the Supplemental Material [24( in more detail. 

For a weak disorder, e* -C T d , the localization length 
at characteristic energy e ~ T d , given by [25H27| 


C(e) 


he 


V/^J 2 ' 


( 6 ) 


greatly exceeds the mean interparticle distance n 1 . This 
ensures the presence of a small parameter: 


-(ft) 


3/2 


< 1 . 


( 7 ) 


In the limit of 7 —^ 00 the fermions are free [17], and 
an arbitrary weak disorder localizes them irrespective of 
their energy, leading to an insulator state at any tem¬ 
perature. A weak interaction between the fermions (and 
thus a strong but finite interaction between the bosons) 
causes many-body delocalization, i.e. destroys the insu¬ 
lator at the critical temperature (see Eq.(42) in Q and 
Supplemental Material (ftjj ) 


T P = 


cs t 


c 


IA In I 


( 8 ) 


Here C ~ 1 is a mo del-dependent numerical constant, 
and Sq = 7 T h 2 kp/mC, is the single-particle level spacing 
on the scale of the localization length £. For a weak 
disorder in the ID case £ « hkpT/m, with r being the 
transport time. Thus we have: 


7 tH 

T 


(9) 


It is known that the interaction between ID fermions 


renormalizes the disorder (see Refs. jl 8 l. |28| ): 

ft = V Q (§) " ; |A| « lj T « Td - 


( 10 ) 


Here 1/tq is the particle-impurity scattering rate. In the 
Born approximation (see Supplemental Material [2J]) 


1 

TO 


2maU 2 
nnh 3 


£* 

h 


2 e* A 

7 ft ) 


1/2 


( 11 ) 


The dimensionless fermionic coupling constant is de¬ 
fined El as: 


A = 


[H(0) - V{2k F )]m 


2'K 2 h 2 n 


( 12 ) 


where V ( q) is the Fourier transform of the interaction 
potential Vp 0 for the transferred momentum q. In the 
secondly quantized form the potential 0 writes: 

= h ^ ^^ki,k 2 ) 

ki,k2,q 

a\ki + q)a){k 2 — q)a(k 2 )a(ki), (13) 

where the operator a{k ) annihilates a fermion with mo¬ 
mentum fc, and 

V{q,k i,fc 2 ) = —- k 2 +2q)(k x - k 2 ). (14) 

m^g 

For q = 2kp the phase-space constraint fixes k\ = — k 2 = 
—kp , and for q = 0 we have (£q — k 2 ) = ±2 kp. This 
leads to [H 


v(2k r ) - n») = SHVf - 8 ' 2,il " 


m 2 g 


m'y 


(15) 
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Substitution of Eq. (fl5l) into Eq. (fl2l) yields the fermionic 
coupling constant 

A = —— < 0. (16) 

7 

Thus, according to Eos. (fTOll and (1S26I) the renormal¬ 
ized impurity backscattering rate and the effective level 
spacing decrease with temperature due to the inter¬ 
fermion interaction 0 - Using Eqs. o, m, and © 
we find: 

Inserting Eq. m into Ea. (IS26l) and using Eq.® we ob¬ 
tain a critical temperature in the quantum degenerate 
regime: 



Equation (flSl) indicates that T c vanishes at 7 = 70 « 8 . 
Note that the weak disorder renormalization group ap¬ 
proach of Giamarchi-Schulz [14] yields the critical Lut- 
tinger parameter K = 3/2 corresponding to 70 = 7.9, 
so that at T = 0 an infinitesimally small disorder leads 
to localization (insulator state) for K < 3/2 (7 > 7.9). 
In this sense, our equation (1181) is in agreement with the 

zero temperature result of Ref. [13- 

In the derivation of Eq. flTSl) we took into account the 
renormalization of the disorder due to interaction of 
the effective spinless fermions with Friedel oscillations 
and neglected the renormalization of the interaction by 
the disorder. It is known H3 that this approximation 
works as long as (7 — 70 ) V 1 / 2 . Close to the zero- 

temperature transition this condition is violated and one 
has to use the coupled renormalization group equations 
for the interaction and disorder. For 7 < 70 this leads 
to the Berezinskii-Kosterlitz-Thouless (BKT) criticality 
with temperature (see the Supplemental Material 0 ) 

T c ~ T d exp[—16?r/ \Jc\V - c %{7 — 70 ) 2 ], (19) 

where ci ~ 8.5 and C 2 ~ 0.93. As soon as 7 be¬ 
comes larger than 70 and (7 — 70 ) exceeds I? 1 / 2 , the low- 
temperature behavior obeys Eg. (TIBI) . 

The derivation of Ea. flTTJl) is similar to that for the 
localization length scale Li oc in the case of spin-1/2 
fermions in Ref. 0. The relation T c ~ Hvf/Li oc , where 
vf is the Fermi velocity, is the same as Eq. m, except for 
the coefficients c\ and C 2 because the critical Luttinger 
parameter for spin-1/2 fermions is equal to 3. However, 
this T c was treated as a temperature of the crossover from 
low to high resistivity. We claim that T c is the insulator- 
fluid phase transition temperature. 

In the limit of extremely strong coupling, 7 > 8 , the 
power of 7 in Eq. m becomes equal to 1 , which gives 

T c ~T d ^. ( 20 ) 

m 7 


The substitution of 7 ~ 1/|A| from Eq.® into Ea. (l20l) 
transforms it to the result of Ref. 3 for T c of weakly 
interacting spinless fermions. Note that in this regime 
the critical temperature does not depend on the sign of 
the interaction between the fermions. 

For a non-degenerate gas where T d -C T, the transition 
temperature can be found by using the single-particle 
picture of Ref. @ after mapping strongly interacting 
bosons onto weakly interacting fermions. We then have: 

5((T C ) ~ nV(T c ). (21) 

The level spacing is S^(e) ~ l/i/(e)£(e), where v{e) ~ 
y/m/h 2 e is the density of states, and the localization 
length is given by Eq.®. According to Ea. (fI41) we have 
V(s) ~ h 2 e/mg , so that nV(T)/T ~ I /7 <C 1. This 
leads to 

T c ~ T d (V"/) 2 / 3 ~ £ » 7 2 / 3 ; T > T d . (22) 

The value of 7 corresponding to a crossover from the 
classical to quantum degenerate regime is obtained from 
equation (1221) setting T c ~ T d . This yields 7 ~ 1/V^> 1. 

It should be noted that at temperatures of the order of 
T c (1221) the fermions remain weakly interacting. The scat¬ 
tering phase shift Ap(/c) in the two-body problem with the 
interaction potential 0 is 

<$ F (fc) - — arctan ( ^ ^ \ _ _ arc ^ an /V (23) 
\2mg/ V 2771 / 

For T < Td the fermions are degenerate and k ~ fcp ~ n. 
The phase shift 5 F (fc) ~ — arctan( 7 _1 ) is thus small as 
long as 7 1, i.e. the fermions interact weakly. In 

the non-degenerate case, T T d , the momentum k is 
of the order of the thermal momentum yjmT/H 2 and at 
very high temperatures the fermions interact strongly, 
Sf ~ 1. However, at T < T c (l22|) the phase shift is still 
small, Sf < 2 ? 1 / 3 / 7 2 / 3 <C 1, and the MBLDT single¬ 
particle picture leading to equation (l 22 l) is valid. 

In Fig.l we show the phase diagram in terms of the 
amplitude of the disorder ( V ) and temperature at a fixed 
interaction strength ( 7 ). As expected, the critical disor¬ 
der for the fluid-insulator transition vanishes as T —> 0 if 
7 > 8 , whereas for 7 < 8 it remains finite. 

The phase diagram in terms of dimensionless temper¬ 
ature T/Td and interaction strength 7 at a fixed weak 
disorder (V <C 1) is displayed in Fig.2. It combines 
the results of the present paper for strong coupling of 
the bosons, 7 > 1 , with the weak coupling (7 1 ) re¬ 

sults of Ref. 3 ■ We clearly see that in a weak disorder 
(T> -C 1) at any finite temperature one has two insulator- 
fluid transitions as 7 increases from small to large val¬ 
ues. The first insulator to fluid transition occurs when 
the interaction between the bosons is weak (7 <C 1). The 
results of Ref. [7] for this case can be written in terms of 
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V 



b) 

FIG. 1. Phase diagram in terms of the disorder strength V 
and temperature at a fixed dimensionless coupling constant 
7 . In a) 7 1. The yellow part of the curve corresponds to 

Eq. (124ll . the blue part to Eg. (12511 . and the red part to Eg. (12611 . 
In b) 7 > 70 ss 8 , and the blue part of the curve corresponds to 
Eo. lfTHll . whereas the red part to Eq. m. Inset: dependence 
of the critical disorder on 7 at T —► 0, described by Eq. HU 
for 7 1 and by Eq. HU for 7 close to 70 « 8 . The dashed 

part of the curve shows the expected behavior for T> ~ 1 . 


the critical coupling 7 c i of the transition as: 

7ci~ D 2 / 3 ; T < V^Ta, (24) 

7 ci~ VT d /T- V x ^T d « T < T d , (25) 
7ci~ V{T d /T )V 2 ; T > T d . (26) 

The reentrance to the insulator takes place at strong in¬ 
teractions (7 O 1 ): 

7c2= 7 o - cVV; T 0, (27) 

7c2 ~ V~ l {T/T d f/ 2 - T»T d . (28) 


For finite temperatures T -C T d the critical coupling is 
determined by a transcendental equation following from 
Ea. ([T 8 l) . The T —> 0 asymptotics corresponds to the 1+1 
BKT quantum transition only slightly modified by the 
weak disorder. In the limit of T> —> 0 equation (1271) thus 
reproduces the result of Giamarchi and Schulz [141 ] . 


T/T d 



FIG. 2. Fluid-insulator transition of finite temperature repul¬ 
sive bosons in a weak disorder, 1 (see text). 


The difference 7 c2 — 7 c i between the two critical in¬ 
teraction constants should decrease as the disorder (V) 
increases. The two insulating regimes at T = 0 are prob¬ 
ably merging ( 7 c i = 7 c2 ) at a certain critical value of the 
disorder, T> c ~ 1 . 


It is feasible to verify experimentally the full picture of 
the finite temperature behavior of ID disordered bosons 
constructed in the present paper. A suitable candidate 
would be the gas of 7 Li atoms where the coupling con¬ 
stant g can be varied by Feshbach resonance from very 
small to very large values 13011 . and the ID regime has 
already been achieved :3lj, The regime of strong 

interactions can be reached, in particular, by using a 
confinement-induced resonance as in the cesium exper¬ 
iments [32| , and the disorder can be introduced by using 
optical speckles like in the first experiments on the ob¬ 
servation of Anderson localization [9j. The insulator-fluid 
transition can be identified in expansion experiments (see 
the discussion in Ref. 0 ), or by analyzing the momen¬ 
tum distribution like in recent experiments for bosons in 
ID quasiperiodic potentials 00]. 


Note Added 

After the present work has been finished, we learned 
about Ref. |34j where the two insulator-fluid transitions 
were observed in the experiment with fermions in the 
one-dimensional quasiperiodic potential. 
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for fruitful discussions and acknowledge support from 
IFRAF and from the Dutch Foundation FOM. The re¬ 
search leading to these results has received funding from 
the European Research Council under European Com¬ 
munity’s Seventh Framework Programme (FR7/2007- 
2013 Grant Agreement no.341197). 
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Supplemental Material: Finite-Temperature Fluid-Insulator Transition of Strongly 

Interacting ID Disordered Bosons 

We detail here five points of the paper. We begin with the description of the disorder potential and the length 
and energy scales of the ID localized states in the absence of interaction, then we explain the fermion/boson duality 
in ID in the presence of an external potential. Furthermore we characterize the finite-temperature fluid-insulator 
transition of weakly interacting spinless fermions in the degenerate and dilute regimes. In the fourth part we detail 
the renormalization of the disorder strength due to weak interfermion interaction. Finally we point out the relation 
between our results and the Giamarchi/Schulz renormalization group formulation of the weakly disordered Luttinger 
liquid. 

Disorder potential length and energy scales. We consider a Gaussian random potential U ( x ) with the correlation 
length (7 and the amplitude Uq (see the main text). The energy-dependence of the localization length is known [sil - [S3j : 

((e) « h 2 e/maUQ , (SI) 

under the condition 

C » <7; (S2) 

so that the extent of the wave-function covers a distance over which the random potential varies many times. Let us 
estimate the energy of a low-lying bound state in the random potential. The kinetic energy of the state is 

K ~ h 2 /2m( 2 , (S3) 

and the potential energy is computed as follows. Consider the segment of the length £, which is made of fragments of 
size <7, where the value of the potential is ~ Z7o - The number of such fragments is (/cr. A wavefunction with the extent 
C gives the potential energy contribution of a typical fragment ~ U^a/C,. In the case of the Gaussian fluctuations the 
typical potential energy of the state is the energy contribution of each fragment multiplied by the square root of the 
number of fragments \J C/cr, that is 

U ~ -U 0 a/<; x JcTo = -C/o^TO (S4) 

Minimizing the total energy E = I\ + U with respect to and taking into account Ea. dSlI) one obtains the length 
and energy scales of the low-energy bound state 

c* = (^/mVf/o 2 ) 1 / 3 , e* = (mer 2 f7o /h 2 ) 1 ^ 3 . (S5) 

Inserting into Ea. (IS'2l) directly leads to the condition stated in the main text: 

U 0 < h 2 /ma 2 , (S6) 


and one then clearly sees that e* -C Uq. 

Boson/Fermion correspondence in ID. Here we extend the arguments of Ref. SJ] to the case of an external (random) 
potential. We demonstrate that although the model is not exactly solvable the ID fermion-boson duality remains 
valid. 

Consider the Schrodinger equation for the 7V-fermion wavefunction x( x i> ■ ■ ■ > x n)' 


(Hi + V F )x(xi,. ■ .,x N ) = Ex(x 1 ,.. .,x N ). 


(S7) 


Here H 1 is the single-particle part of the Hamiltonian, which differs from Eq.(3) of the main text by the presence of 
the external (disorder) potential U(x): 


Hi 




(S8) 


The two-body interaction part of the Hamiltonian, Vp. is determined by Eq.(4) of the main text, which can be 
rewritten as 

A fe4 

V F = — 2 - ^2 djk$(xjk)d jk , (S9) 

^ l<j<k<N 
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where x jk = Xj - x k and d jk = (d Xj - d Xk )/2. 

For non-coinciding coordinates Eq. rtSTl) reduces to 

■ -,xn) = Ex(x i,.. .,x N )- (S10) 


This equation should be supplemented by the boundary conditions, which determine the behavior of the wavefunction 
in the limits x 3 —> x k ±{). In order to find these conditions one integrates Ea. (IS7l) over an infinitesimally small interval 
—x < Xj k < x; x, x —> +0 of Xj k . Using Eas. dS8l) and (IS9l) and taking into account the identity d 2 . +d Xk = 2d 2 k + \d\. k 
(here Xj k = (xj + x k )/2) we then obtain: 



4ft 4 

_l - 2~ 

m z g 


dj k S(x jk )d jk \x{x jk ,...) + Regular = 0, 


(Sll) 


where “Regular” includes terms that vanish in the limit x, x —> 0. Performing the integration in the first term we 
obtain the condition which does not contain the external (random) potential: 

h 2 Ah' 1 

- d x [x(x, ■■■) - x(~x, ■ ■ ■)] d- 2 ~[S(x)d x x{x,...) - 5(-x)d x x(~x , ...)]= 0. (S12) 

m m z g 

In the limit x,x —> +0 equation (IS 121) can be written as 


kX(.Xj k +0, • • •) d jk x{x jk 0,...). (S13) 

Integrating Ea. (IS12[) over x in an infinitesimally small interval —x < x < x and taking the limit x — > +0 we have: 

Ah 2 

X(xjk = +0,...) - x(xjk = -0,...) =- djkx{xj k , ■ ■ • )|^ fc =+o- (S14) 

Trig 

The fermionic wavefunction is antisymmetric with respect to the transmutation of Xj and x k . Accordingly, xi x jk = 
+0,...) = — x(xj k = —0,...), and in terms of Xj,x k equation (IS14() takes the form: 

(d Xj - d Xk - ■ - ,Xj,Xk,...)\ Xj =x k + o = 0. (S15) 

Note that Eas. dSlOl) and (| S1 -5 [) could equally be obtained for the Lieb-Liniger Hamiltonian of Eq.(l) of the main 
text, completed by the external potential: 


H b = Hi +g^26(xj - x k ), (S16) 

j<k 

which acts on the N-body bosonic wavefunctions. As a consequence the many-body spectra of the models (1571) 
and (IS16I) are identical and the eigenfunctions are equal in magnitude and only differ by a sign upon exchange of 
coordinates. 

Conductor-insulator transition of disordered weakly interacting ID spinless fermions with gradient-dependent inter¬ 
action. In the case where all single-particle eigenstates are localized (as in ID), there is a critical temperature T c 
above which weakly interacting spinless fermions turn to a conducting state. Here we characterize the transition both 
in the quantum degenerate regime jsd| and in the high-temperature dilute regime [S6j. In the basis of the localized 
single-particle states ip a (x) with energies s a the many-body Hamiltonian Hi + Vp , with Vf given by Ea. (IS9l) . writes 

Hf = i ^2 VfffpO^alapan, (S17) 

ot otfi'yS 


where a a are the field operators and the two-body interaction matrix elements Vffp are given by: 


V 1& — - 

V ap ~ 


The single-particle states take the form 


4ft 4 /' 

2 / d xif*{x)[d x tpg{x)][d x ifp(x)]'if a (x). 

m z g J 


-Ix-Xcj/Cc 


(S18) 


ipa(x) = A a (x)e 


(S19) 
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where A a ( x) is an oscillating function with wave number 

k a = \/2m£ a /h. 


Therefore the matrix elements (IS18I) evaluate to 

v 2p ~ [ dxA 7 (x)A s {x)A 0 (x)A a (x)t 

m 9 J 


— (\x—Xri\/C,~t + \x—X5\/C5 + \x-Xp\/Cp-\-\x—X a \/Cat) 


Taking into account the condition that the levels a ,7 and /3,<5 are neighbors in energy [H,@] , we obtain 

48 H 2 £ B 


V 


7<5 


m 2 g max(Ca,C/3) mg max(C«, C/ 3 )' 


(S20) 


In the quantum degenerate regime, T -C Td = n 2 h 2 n 2 /2m (n is the average density of particles), the two-body 
processes responsible for the fluid-insulator transition occur at the energy scale Td- The relevant two-body matrix 
elements are 


Kl ~ Ad c , (S21) 

where the effective interfermion coupling parameter A is of the order of the inverse of the Lieb-Liniger parameter: 

A~ —e-«1, (S22) 


h 2 n 1 

-= - < 1 , 

mg 7 


and the spacing between the single-particle energy levels in a segment of length £ writes 


(S23) 


with v = m/ir 2 h 2 n being the density of states at the Fermi energy. The many-body Hamiltonian (IS 1 71) with the 
two-body matrix elements (IS21I) leads to the critical temperature in the quantum degenerate regime [S5| : 


T r = 


CSt 


|Aln|A||’ 


(S24) 


where C is a mo del-dependent numerical constant of the order one. Taking into account the duality (IS22I) and the 
level spacing renormalization due to interparticle interaction (described in the next part), we obtain Eq.(18) of the 
main text. 

In the dilute high-temperature case T Td, the average occupation of the single-particle energy levels is small and 
given by the Boltzmann distribution. In that case the temperature T gives the relevant energy scale of the states 
participating in the interaction processes. The critical temperature in this regime is given by [S 6 ( 


with the one-particle level spacing at energy e: 


d c (T c ) ~ nV{T c ), 


6 c {£) = 


zz(e)C(e)' 

The energy-dependence of the density of states is 

v{e) = 1 ym/2n 2 H 2 £, 


(S25) 


(S26) 


(S27) 


and the localization length £(e) is given by Ea. dSlI) . In Ea. (|S25l) the quantity V(e) is the Fourier transform of the 
two-body interaction potential at momentum k = y/2m£/H: 


V(e) 


4 H 2 £ 
mg ' 


(S28) 


Inserting Eas. dSTI) . (IS27I) . (IS26I) . and (IS28I) into (IS25I) we obtain Eq.(22) of the main text. 
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Renormalization of the impurity scattering amplitude. In this section we describe the renormalization of the disorder 
strength by the interaction and derive the resulting Born approximation backscattering rate in the case of many 
impurities. The renormalization of a single impurity transmission amplitude t due to weak interaction between 
spinless fermions was computed by Matveev et al. |S7j |. The result is a temperature dependence of the transmission 
coefficient T = \t\ 2 : 


T(T) = 


1 


(S29) 


l + (l/T 0 -l)(T d /T) 2 *’ 

here 7o is the impurity transmission coefficient in the absence of interactions. The coupling parameter is defined as 


y(o) - v(2k F ) 


(S30) 


with V{k ) being the Fourier transform of the interfermion potential, and vf = irhn/m the Fermi velocity. 

In quantum mechanics |S8| a particle of energy T d has an amplitude of transmission through a potential of extent 
a and amplitude Uq given by 


1/T — 1 = mo 2 U$/2H 2 T d 


(S31) 


under the condition (IS6I) . On the other hand the backscattering rate caused by many impurities with concentration 
Cimp is evaluated in the Born approximation (Tt h) as 


and taking into account Ea. (lS31l) we get 


1/r = 2c irnp (j 2 UQ/h 2 v F: 


l/ T= 2nhnc ' m P (l/T-l). 


(S32) 


(S33) 


The prefactor in the right-hand side of Ea. (IS33p does not depend on the properties of the single impurity potential 
and therefore it is not affected by the renormalization (IS29I) . As a consequence we can write the renormalized 
backscattering rate as 


1/r = (1 /to)(VT) 


2A 


(S34) 


where setting c, mp = 1/tr the bare backscattering time is 


nh 3 n 

2mcrUQ 


(S35) 


Giamarchi and Schulz renormalization group for the disordered bosons. We now turn to the Luttinger liquid 
formulation of the disordered ID bosons and point out the relation between its predictions and our theory. The 
renormalization group approach of the disordered Luttinger liquid of bosons was developed by Giamarchi and Schulz 
m an d it describes the zero-temperature algebraic-superfluid to Bose-glass transition at mo derat e interactions and 
small disorder. The influence of disorder is taken into account by the backscattering field £(:r) IH3 , with correlations 

(x')) = TT 2 h 2 v 2 nT>S(x — x'), (S36) 


where the symbol (.) represents the average over disorder realizations, v = vf/K is the velocity of the sound mode, 
K is the Luttinger parameter, and T> <C 1 is the dimensionless parameter of the disorder. The relation between V> 
and the parameter T> that we use in the main text is 



7r 



K 2 
2 3 / 2 7T 


V. 


(S37) 


The scaling transformation of the disordered Luttinger liquid model leads to two coupled differential equations 
giving the flow of K and T> with respect to the scale parameter. The solution of the equations leads to the emergence 
in the localized phase of the length scale Li oc which diverges at the glass-superfluid transition. There are two limiting 
cases where Li oc takes on two different forms. When K < 3/2 and 2? 1 / 2 -C 3/2 — K this length scale writes 

L loc ~a( i/p)37abr, (S38) 
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where a is a short-distance cut-off that we set equal ton 1 . Thus for infinitesimally small disorder the theory predicts 
the fluid-insulator transition at K = 3/2. At this point the parameter of the disorder strength is 

V = « 0.252?. (S39) 

Having obtained the length scale Li oc one can infer the critical temperature mu 

T c ~ hv F /L loc . (S40) 

It is interesting to note the correspondence between Ea. (lS40|) and our strong coupling critical temperature in the 
degenerate regime given in the main text. Indeed given the development of the Luttinger parameter with respect to 
l/q = h 2 n/mg -C 1: 


K = 1 + 4/7 + 0(l/7 2 ), (S41) 

we see that our case corresponds to the first-order strong coupling expansion of K and our result extends Ea. (lS40|) 
to the whole temperature regime T < Td = Tr 2 h 2 n 2 /2m. The two results are in agreement in the regime 2 ? 1 / 2 -C 
(7 ~ 7 o)/ 7 o < 1 an d the upper limit is equivalent to T/Td <C 2?, where 70 « 7.91(1) is the solution of the equation 
AT ( 7 ) =3/2 evaluated numerically from the Bethe Ansatz equations mi- 

On the other hand when 2 ? 1//2 >> K — 3/2 ( K >3/2) the renormalization of the interaction is important and we 
must solve the coupled differential equations for K and V. This leads to the length scale 


Li oc ~ a exp 


7r/y 92?/8 — (K — 3/2 ) 2 


(S42) 


and accounting for Ea. (IS40l) and the dependence K( 7 ) near K = 3/2: K( 7 ) « 3/2 + 0.058(7o — 7 ), we obtain Eq.(19) 
of the main text. 
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